■ DIFFUSION IN ONE DIMENSIONAL RANDOM MEDIUM 

Q\ 

AND HYPERBOLIC BROWNIAN MOTION 

C ' 
i— » 

^° . 

Alain Comtet 1 and Cecile Monthus 2 

(N 
O 

■ Division de Physique Theorique 3 , IP N Bat. 100, 91406 ORSAY, France 

O 

in 

■ and 

o 
o 

> 

•i-H . 

' Abstract 

&. 

Classical diffusion in a random medium involves an exponential functional of Brownian 
motion. This functional also appears in the study of Brownian diffusion on a Riemann 
surface of constant negative curvature. We analyse in detail this relationship and study 
various distributions using stochastic calculus and functional integration. 

Keywords : Brownian motion, random media, negative curvature. 



IPNO/TH 95-18 

1 comtet@ipncls.in2p3.fr 

2 monthus@ipncls.in2p3.fr 

3 Unite de recherche des Universites Paris 11 et Paris 6 associee au CNRS 



1 



1. Introduction 



There is a close link between one dimensional, or quasi-one dimensional, disordered 
systems and Brownian diffusion on Riemann manifolds of constant negative curvature. 
Such a correspondance can be traced back to the pioneering work of Gertsenshtein et al. 
W who have shown that statistical properties of reflection and transmission coefficients 
of waveguides with random inhomogeneities are directly related to some random walk 
on the Lobachevsky plane. There has been a renewed interest in this approach through 
the study of mesoscopic systems. The description of quasi-one dimensional mesoscopic 
wires involves a Fokker-Planck equation giving the probability distribution of the N 
eigenvalues of the transmission matrix ( 2 ) . Recently it has been shown that this equation 
can be interpreted as the diffusion equation on a Riemannian symmetric space ^ . 

The purpose of this work is to show how the one dimensional classical diffusion of a 
particle in a quenched random potential U (x) that is itself a Brownian motion, possibly 
with some constant drift, is directly related to Brownian motion on the hyperbolic plane. 
Since the latter is the archetype of chaotic systems ^ , our work forms a bridge between 
disordered and chaotic systems. 



2. Fundamental random variable for one dimensional classical diffusion 
in a quenched random potential U(x) 

A large amount of work has been devoted to random random walks, defined on a 
lattice by the following Master equation 

Pn(t+1) (t) + (3 n+1 P n+1 (2.1) 

where a n is the random quenched transition rate from site n to site (n + 1), and 
(3 n = 1 — a n is the random quenched transition rate from site n to site (n — 1). From 
a physical point of view it is convenient to introduce a corresponding random potential 
U (n) on each site n and to write the ratio of the two transition rates a n and f3 n as an 
Arrhenius factor 

a /3[u(n-l)-U(n)] 

Pn _ g _ J3[U(n+l)-U( n -l)] (2 2) 

The study of different physical quantities related to this random random walk ( 5 ) involves 
systematically random variables of the following form 

b n 

Z(a, b) = ^ II a k = °a + Va°a+l H V Va°a+l ■ ■ ■ °b (2.3) 

n=a k=a 
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Their fundamental property is to satisfy the linear random coefficient recurrence relation 



Z{a, b) = a a 



l + Z{a+l,b) 



(2.4) 



This discrete multiplicative stochastic process also appears for the Ising chain in a 
random magnetic field ( 6 -*. 

Let us now consider the continuous model of classical diffusion on the line defined 
by the Fokker-Planck equation for the probability density P(x, t \xo, 0) 

dP 1 d fdP 



where {F(x)} is a quenched random force. In this continuous limit the discrete random 
variable Z(a, b) defined in (2.3) becomes the following exponential functional of the 
random potential U(x) = — f x F(y)dy 

r(a, 6) = f to F(»> " "(»)] = f * f ^ (2.6) 

J a J a 

The evolution of this functional is governed by the stochastic differential equation 

^-(a,b)=PF(a)r(a,b)-l (2.7) 

which replaces the random coefficient recurrence relation (2.4) satisfied by Z(a, b). Note 
that the stochastic term F(a) appears multiplicatively, so that the fluctuations of the 
random force are coupled to the values taken by the random process r(a, b). 

Let us now make clear how the functional r(a, b) arises in some physical quantities 
associated with the classical diffusion of a particle in the quenched random environment. 

• If the random force has a positive mean < F(x) >= Fq > 0, the probability distribu- 
tion of the random functional 

= lim r(a,b) (2.8) 

(b-a)^oo 

determines the large time anomalous behavior of the Brownian particle position ( 7 ) . In 
particular, the velocity defined for each sample as 

d f +c ° 

V= lim — / dx xP{x,t |x ,0) (2.9) 

t^oo dt J_ 00 

is a self-averaging quantity inversely proportional to the first moment of Too ^ 



When the quenched random force is distributed with the Gaussian measure 



~ J dx [F(x)-F ] 
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VF(x) e Za J (2.11) 
the probability distribution Voo(t) of the functional reads ( 7 ) 

Poo(t) = =7 — r- ( — ) e ctr oc — — (2.12) 

where = ^ L > is a dimensionless parameter and a = The algebraic decay at 
large r accounts for the dynamical phase transitions which occur in this model ( 7 ^. In 
particular, equation (2.10) implies that the value \i = 1 separates a phase of vanishing 
velocity V = for < /j, < 1, and a phase of finite velocity V > for ji > 1. 

• The functional r(a,b) also arises in the study of transport properties of finite size 
disordered samples. The stationary current Jn which goes through a disordered sample 
of length N with fixed concentrations Pq an d Pn & t the boundary can be written in 
terms of the exponential functional tn = r(0, N) as ( 8 ) ( 9 ) 

1 



Jn= 2 



P „ d\lYT N 
— ±N 



(2.13) 



r N JV 

When the end x — N is a trap described by the boundary condition P/v = 0, the flux 
is simply a random variable inversely proportional to tjv- The probability distribution 
of tat has been studied for the case of zero mean force Fq = and for the general 
case with arbitrary mean force ( 9 ) by different methods. 

It is also interesting to point out that a lot of mathematical work has recently 
been devoted to the functional tn in relation with finance ( 10 - ) ( n -* ( 12 -*. 
In ( n ), Yor pointed out the relation between the functional tn for the particular case 
fj, = 2> and free Brownian motion on the hyperbolic plane. In the following, we first 
rederive this correspondance, and then generalize it to arbitrary fj,, using some external 
drift on the hyperbolic plane. 



3. Relation with hyperbolic Brownian motion 



The upper half-plane { (x,y) , y > 0} endowed with the metric 

ds , = dj^+ds 

defines a two-dimensional Riemann manifold of constant negative Gaussian curvature 
R = — 1. The surface element dS and the Laplace operator A are covariantly defined as 

dx dy , 9 / d 2 d 2 \ , . 

dS = — ? r L ^d A = y 2 — - + _ (3.2) 



y 2 \dx 2 dy 2 



Free Brownian motion on this manifold is denned by the diffusion equation for the 
Green's function G t (x, y) 



d 2 



dG 2 f d 2 

dt ^ \dx 2 ' dy 2 



It is convenient to choose the initial condition at the point (x = 0, y = 1) 

G t (x,y) t -j + 8(x) %-l) 

The normalization of the Green function Gt(x, y) then reads for any time t 

/r+oo r+oo ^ 

dS G t (x,y) = J dx J dy G t (x,y) 



(3.3) 



(3.4) 



(3.5) 



Consider the probability density Pt(x,y) 



1 



Pt(x,y) = -J G t (x,y) 

y 

normalized with respect to the flat measure dx dy 

/+oo r+oo 
dx / dy P t (x,y) 
-oo JO 

This probability density satisfies the Fokker-Planck equation 



dt V dx 2 dy 2 



if d 2 , , n 

+ — I y 2 P 



(3.6) 



(3.7) 



(3.8) 



and the initial condition at t — 



P t (x,y) t -j + 8(x)8(y-l) 



(3.9) 



We now introduce two independant Gaussian white noises rji(t) and 772 (t) and write 
the stochastic differential equations for the process {x(t),y(t)} corresponding to the 
Fokker-Planck equation (3.8) following respectively the Ito or Stratonovich convention 

(13) 

ml at \ (3.10) 

% = V2Dy V2 (t) J 



Stratonovich 



/2Dy m (t) 
-Dy + V2Dym(t) J 



(3.11) 
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Integration of these two coupled differential equations gives 

-Dv + V2D [ r] 2 (s)ds 



x(t) = \p2~D [ dv 771 (u) e 
Jo 

-Dt + V2D / 7] 2 (s)ds 
I y(t) = e Jo 



(3.12) 



The process y(t) is therefore simply the exponential of some linear Brownian motion 
with negative drift. Let us now study more precisely the process {x(t)}. Since the white 
noise rji appears linearly in x(t), all the odd moments therefore vanish 



< x 2n+1 (t) > m =0 



(3.13) 



and all the even moments can be rewritten as averages over r]2 only. In particular the 
second moment reads 



<x 2 (t) > VuV2 =2D < 



-2Dv + 2V2D 



dv e 



Jo 



s)ds 



> 



If we set 







[ F(u)du = 
Jo 



2Dv-2V2D / rfr(s)ds 



(3.14) 



(3.15) 



the expression between brackets on the right hand side of Eq (3.14) is nothing but the 
functional 

-P [ F{u)du 



n 



dv e 



(3.16) 



which is encountered in the study of classical diffusion in a quenched force {F} dis- 
tributed as a Gaussian white noise 



V F(x) e 



1 r + °° 



'dx 



(3.17) 



with parameters F = and a = 8 More generally, all even moments of x(t) are 
proportional to moments of the functional Tt 



< x 2n (t) > m , m = (2Dr < t? > F M{n) 



(3.18) 



J\f{n) being the combinatorial factor coming from Wick's theorem that counts the num- 
ber of ways to pair the 2n fonctions rji. In particular, Af(n) is equal to the moment of 
order 2n of a suitable Gaussian random variable £ with variance unity 



N{n) = 



+ °° dP 

^ i 2n e 2 = 2 



2tt 



O _on r ( n til = -Q (2 fc_l) 



(3.19) 



k=i 



These identities between moments allow us to rewrite the process x(t) in the very simple 
form suggested by Yor ( n ) 

x(t) = V2D / du rj(u) (3.20) 
Jo 

where rj is a new independent Gaussian white noise. The process x(t) is therefore 
identical in law to a linear Brownian motion whose effective time is the random variable 
Tt- This result can also be obtained through a functional integration method using an 
appropriate change of time in path-integral (see Appendix). 

To sum up, the free Brownian motion {x(t), y(t)} on the hyperbolic plane can be 
rewritten in terms of the two independant white noises of measure 

-tt: [F(x)-F fdx -L dt V \t) 

V F(x) e 2a Vrj(t) e 2 J (3.21) 



The process y(t) is simply the exponential of the Brownian motion with drift 
U(t) = -f*F(u)du 

y(t)= e 2 Jo = e2 [ ) (3.22) 

The process x(t) can be viewed as a linear Brownian motion 

x(t) = V2D / du 7](u) (3.23) 
Jo 

with an effective time r t that is itself a random process depending on y(t) 

T t = f dv eP U ^ = f dv y 2 (v) (3.24) 
Jo Jo 

Note that the dimensionless parameter \i = ^0 which characterizes the different phases 
of anomalous diffusion is \ for the free hyperbolic Brownian motion. This value rep- 
resents the natural drift induced by the curvature of the Poincare half-plane. It is, 
however, easy to generalize the previous analysis to arbitrary fx with the introduction 
of some external constant drift m in direction y 

/j=^+m (3.25) 
The corresponding stochastic differential equations then read 

^ = V2Dy Vl (t) 
It6{ } (3.26) 

^ = -2Dmy + V2Dy m (t) 
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Stratonovich 



dx i ~ . . 



dy 



dt 



-2Dvy+V2Dy m (t) j 



(3.27) 



For any \i 1 there is therefore a direct correspondance through Eqs (3.22-3.23) be- 
tween the joint stochastic process characterizing the one dimensional diffusion { random 
potential U{t), exponential functional r t } and the Brownian motion {x(t), y(t)} on the 
hyperbolic plane with possibly some external constant drift m along direction y. We 
now consider some consequences of this correspondance. 



4. Marginal laws of the processes r t , x(t) and y(t) 



The marginal law Y t (y) of the process y(t) reads according to eq (3.22) 

Y t (y)= VU(s) e 2a Jo \ds J 6 I y - a e% W (4.1) 

Ju(o)=o y J 

We therefore get after some algebra the following log-normal distribution 

1 2 
1 -t^-TWj/) + 2uDt] 

Y ' w ^ wsm e iDt (4 - 2) 

where \i = ^0- and D = In the case of free Brownian motion [ji = |), this marginal 
law tends to a 5 distribution in the limit t — > oo 



Y oc (y) = S(y) (4.3) 

The Brownian particle is therefore attracted to the y = result of the curva- 

ture of the hyperbolic plane. Note that this axis represents infinity on this plane. This 
limit law remains unchanged as long as fx = [\ + m) > 0. However, when the constant 
external drift m along y is negative enough to overcome the natural drift of the hyper- 
bolic plane (m < — |), there is no equilibrium distribution for the process y(t). In fact 
for n > 0, the existence of this stationary distribution for the process y(t) governs the 
existence of a stationary distribution for the process x(t) that we now construct. 

The process x(t) is a Brownian motion of effective time r t which is a functional of 
the process y(t) 

x(t) = V2D / du 7](u) with r t = / dv e P U ( v ) = dv y 2 (v) (4.4) 
Jo Jo Jo 
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The statistical independence of rj(t) and r t allows us to write the marginal law X t (x) 
of the process x(t) in terms of the probability distribution V , t( T ) of the functional r t 



X 

e~4Dr 



(4.5) 



In a previous work ( 9 \ we showed that V , t( r ) satisfies a Fokker-Planck equation, associ- 
ated with the stochastic differential equation (2.7) which in the Stratonovich prescription 
reads 

dMr) _ d 



dt 



dr 



«r 2 ^^+((^ + l)ar-l)^(r) 



(4.6) 



It has to be supplemented by the initial condition ip t= o(r) = S(r). We recall the solution 
found in terms of an expansion in the Fokker-Planck eigenvectors basis 



0<n<f 



T(l + /j-n) \ar J n \ar J 



> tn( »- n) (-inp-2n) ( 1 \ 



l + M 



4tt< 



Jo V 2 2/ \ar ) \ar J 

(4.7) 

This leads after integration with the Gaussian kernel (4.5) and the transposition a = AD 
to an expansion presenting the same time relaxation spectrum 



Y(<A- V r - ADn ^ ~ n )t !~12^ ~2n)r(^+i-n) / 1 \ 
t[) Jt* n!r(l + /i-n) r(i-n) \l + x* J 



0<n<f 



x F ( —n, n- ix, -; -a;' 



,V Has e-^+^ssinhns r (-£ + 
47r 2 v^r 7 V 2 



I h z- 

2 2 



+ 1 S jtt + 1 .S 1 2 

x F h z-, %-, -; -x 2 

V 2 2' 2 2' 2' 

where F(a, 6, c; z) denotes the hypergeometric function of parameters (a, 6, c). 
For /i > 0, there exists an equilibrium distribution X 00 (x) 
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(4.8) 



(4.9) 



The existence of this limit law for hyperbolic Brownian motion reveals a "localization" 
phenomenon in direction x. This effect comes from the attraction towards the axis 
y = 0. Note that for the free case {ji = \) the asymptotic marginal law X 00 (x) is 
simply a Lorentzian 

x - (x) = lih (4 - 10) 

In summary, there exist equilibrium distributions for and -X"oo 

as long as 

H = {\ +m) > 



Y 00 (y) = 6(y) 



i r{fjL + i) 



r(/x) \i + x 



} (4-11) 



The joint law of the processes {x(t),y(t)} therefore presents the following factorized 
form in the limit t — > oo 

P 00 (x,y) = X 00 (x)S(y) (4.12) 

However, as long as time t is finite, the two processes x(t) and y(t) remain coupled. 
The study of their joint law is then needed to get a complete description of hyperbolic 
Brownian motion. 



5. Joint laws of the processes {r t ,U(t)} and {x(t),y(t)} 

Path-integral methods allow us to obtain very simply the joint law iptij \\ u) of 
the random variables 

ft 

r t = dx e W x ) and U(t) 
Jo 

Let us generalize our previous approach ( 9 ) . The r-Laplace transform of the joint law 

E (e~P Tt W u )=J o dr e ~ PT Mr II u) (5.1) 
can be written as a path-integral over the random potential 

, \ f u(t)=u -1- f (^L + fX dx-pfdxeWW 

E (e-P Tt \\u)= VU(x) e 2(7 Jo \dx J J 

F °* F ° u /•"«>=« -— — \ J*-« I' dxrM') 



2a Jo \dx J Jo 



2a e a / VU{x) e la Jo \ ax J j ( 5-2 ) 

Ju(0)=0 
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The remaining path-integral is simply the Euclidean quantum mechanics Green's func- 
tion < -u|e~* H |0 > associated with the Liouville Hamiltonian 

We therefore get 

Fit F 

£ ( e~P Tt || = e 2a e a <u\e~ tH \Q> (5.4) 

The expansion of the Green's function < -u|e~*^|0 > in the basis of eigenf unctions ipk(u) 
of the Hamiltonian H 

^ 2 \/- sinh ^ X2 '^ ( 2 ^ e "" /2 ) (5 ' 6) 

gives 

/■+°° AU „ 

< w |e- tH |0 >= y — ^(0) e" fc * (5.6) 

We finally obtain that the r-Laplace transform of the joint law (f> t (r \\ y) of the random 
variables 

n=fdxeP U W and y(t) = e2 U{t) 
Jo 

reads after some changes of variables 

E(e-™\\ y) = dr e~P T <j> t {r \\ y) 

at o 

— -u at 9 , . 

e 4 1 ~~r$ -i r . / /p 



7T 2 yi + A 4 



y dge 4 9 ? sinh7r ? K iq (^y^j K iq (^2^j (5.7) 



Let us now compute the x 2 -Laplace transform of the joint law Qt(x,y) of the 
hyperbolic Brownian motion starting from the series of moments of x 2 (t), with y(t) 
being fixed as y 

E (e-** 1 ® || y) = y°° efe Q t (z,y) 

= E ^ E ( x2n M II */(*) = y) ( 5 - 8 ) 

n=0 

Eq (3.23) allows us to write the following relation between moments of x 2 (t) and mo- 
ments of r t when y(t) is fixed as y 

E (x 2n (t) || y) = (2DT M(n)E (r?(t) \\ y) (5.9) 
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_e 

N(n) = I i 2n e 2 (5.10) 



with the combinatorial factor N{ri) introduced in (3.18-3.19). Using the integral repre- 
sentation 

j-oo v/2~7r 

we can resum the series of moments of r t under the integral 

E ( e -9* 2 (*) || y ) = | + °° J= e "T £ (V^£ 2 T t || y) (5.11) 

Using Eq (5.7) and the correspondance a = 4Z) we get 

E ( e -qx\t) ||^ = (512) 



/c 2 

1 / 1 \ ^ 



J dke 4 1 J dve- Dt W +u hsinhw K iv {ky) K iv {k) 



K 2 ^/Tfq \yj 

For the free case (// = ^), this result can readily be recovered from the expression of 
the Green's function G t (x,y) on the hyperbolic plane ( M ) through 



E 



fe-^W || y] = f + °°dx e-^ 2 -1 G t (x,y) (5.13) 

v / j—oo y 



Let us finally mention that an alternative expression of the joint law t (r || y) given in 
Eq (5.7) has been obtained by Yor ( 12 ^ through the time Laplace transform 

r dt e~ st Mr II V) = 7+77 e"(^) /„ (—) where z/ = A A/ 2 + 4- (5.14) 

More generally, we refer to the mathematical literature ^ 10 - ) < - 12 - ) for different ex- 
pressions related with probability distribution of the functional t{T s } where T s is an 
independent time, exponentially distributed with parameter s. 



6. Conformal mapping from Poincare half-plane to unit-disk 

We introduce the conformal mapping from the Poincare half upper-plane {z = 
x + iy , y > 0} to the unit disk {w = re* 9 , \r\ < 1} 

iZ+1 

w — 6.1 

z + i 



12 



The radial coordinate r is directly related to the hyperbolic distance d on the Poincare 
half-plane between the arbitrary point (x,y) and the point {x = 0,y = 1} that we 
choose as the initial point of Brownian motion (3.9) 



r = tanh 



(6.2) 



The circle at infinity r = 1 corresponds to the axis y = 0. The unit disk is very well 
suited to study the free hyperbolic Brownian motion since it contains explicitly the 
rotational invariance in the angle 9. Let us write in the new coordinates the metric 



ds 2 = 

and the Laplace operator 

d 2 



dx 2 + dy 2 

y 2 



(1 -r 2 ) 2 



(dr 2 + r 2 d0 2 ) 



(6.3) 



A = y 2 



+ 



d 2 



dx 2 dy 



2\2 



(1 - r 2 ) 



ld_ fd_\ 1 >i : 



r dr V dr 



36 2 



(6.4) 



In the free case, the Fokker-Planck equation for the probability density Q t {r,9) 



reads 



dQ 
dt 



d_ 

dr 



d f (1 — r 
dr 



2\2 



4r 



-Q 



(l-r 2 ) 2 d 2 

+ 4r 2 QQ2^ 



(6.5) 



The rules of stochastic calculus ( 13 ) give the corresponding stochastic differential equa- 
tions „ „ 

(l-r 2 ) 2 



Ito 



dr =D . 
dt 



4r 



dt 



'2D 



2r 



+ V2D 



Ve(t) 



1 



Vr(t) 



(6.6) 



Stratonovich 



dr 
dt 



1 — r 4 / 1 — r 2 

D-^+VtD—j-Vrtf) 



d9 , l-r 2 , , 



(6.7) 



These equations cannot be integrated straightforwardly to give the processes {r(t), 9(t)} 
as functionals of the white noises {rj r (t), rjg(t)} , unlike the representation (3.12) for the 
processes {x(t), y(t)}. However, the asymptotic probability distribution Qoo(r, 9) in the 
limit t — > oo can be written directly from symmetry considerations. As expected, it is 
simply the uniform measure on the unit circle 



Q 00 (r,9) = ^S(r-l) 



(6.8) 



For n 7^ |, the external constant drift m = (/i — |) along direction y breaks the 
rotational invariance in the angle 9, and the unit disk is not particularly well suited any- 
more. Nevertheless, it is easy to write the asymptotic probability distribution Qoo(r, 9) 
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from the asymptotic law (4.12) for P 00 (x,y) and from the expression of the Jacobian 
J{(r,0)/(x,y)} 

n , m = 4r / 2rcos6> 1 - r 2 \ 

Voolr ' J (r 2 + l-2rsin#) 2 00 Vr 2 + 1 - 2r sin6>' r 2 + 1 - 2r sinfly' 1 J 

We obtain after some algebra the generalisation of (6.8) 

QUr.9) = Hr - 1) ^-^f 1 (1 - sinf)"" 4 (6.10) 



7. Conclusion 



We have explored in great details the relations existing between, on the one hand, 
the exponential functional r t that governs most of the transport properties of one dimen- 
sional classical diffusion in a random potential U (x) , distributed as a Brownian motion 
possibly with some drift, and on the other, free or biased hyperbolic Brownian motion. 
As discussed in the introduction, this has to be considered as a special example of the 
more general link that connects one dimensional, or quasi-one dimensional, disordered 
systems, which usually admit some multiplicative stochastic stucture, and Brownian 
diffusion on symmetric spaces. 

It is interesting to point out that the time relaxation spectrum found for the 
probability distributions ip t {r) and X t (x) (Eqs 4.7 and 4.8) also appears for the quantum 
spectrum of a particle subject to a constant magnetic field B on the hyperbolic plane 
( 15 ). For the hyperbolic geometry, constant magnetic field is defined as the flux per 

covariant surface element dS = —^dxdy. The two spectra coincide if we choose the 

y 

following correspondance between magnetic field B > and drift \i 

2 

In this context, the existence of bound states in a strong enough magnetic field B > | 
corresponds to the presence of closed classical orbits ( 15 ^. It would be particularly 
interesting to understand more deeply this correspondance between spectra at the level 
of stochastic processes themselves. 

As a final remark, let us mention that the Liouville Hamiltonian that we encoun- 
tered in the path-integral formalism (5.3), and that is closely related to hyperbolic 
geometry, also appears in the study of refined properties of one dimensional quantum 
localisation for the Schrodinger Hamiltonian H = — -£? + V(x) where V(x) is a Gaussian 
white noise potential ^ 16 ^. Kolokolov uses a path-integral method to compute correla- 
tion functions of eigenstates and distribution function of inverse participation ratio in 
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the high energy limit. In this formalism, Liouville Hamiltonian shows up in an effective 
action of path-integral. Expansion of this path-integral in a basis of eigenstates then 
gives expressions very similar to the one we obtained for the probability distribution of 
the joint law of the processes {r t , U(t)}. 
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Appendix : Path-integral method to prove the identity in law (3.20) 

In paragraph 3 we derived the identity in law (3.20) mentioned by Yor for the free 
hyperbolic Brownian motion ( n ) through a moments calculation. As explained before, 
this identity can easily be generalized to any fx by the same method. Let us now derive 
it through a more straightforward path-integral method using an appropriate stochastic 
reparametrisation of time in path-integral. This tool has already proven to be very 
useful in other contexts ( 17 K 

The integration of stochastic differential equations (3.27) gives the process x(t) as 
the following functional of the two white noises {rji, 772} 





The marginal law X t (x) of the process x(t) therefore reads by definition 






Let us first change from r)2{u) to U{u) defined by ^U(v) = —2D/iv + y / 2~D r]2(s)ds 



r -- f duril(u) r [ du (^-^ + 2D^\ 

= V Vl (u)e 2 Jo / VU{u)e AD Jo \ 2 dv J 

J Ju(0)=0 
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-U(v] 



S\x-V2DJ dvrn(v)e2 v/ ) (A3) 
And now from 771 (it) to x(u) = \[2D J™ dv rji(v) e^ u ^ 



X t (x) = 



r x(t)=x f -— / du -— + 2Dfi) 

/ Vx{u) / VU{u) e 4D Jo \2dv J 

Jx(0)=0 Ju(0)=0 



'£c(0)=0 J[/(0)=0 

2 



aw 



(A4) 



Let us now perform a time-reparametrisation of the trajectories a; (it) in order to recover 
the Wiener measure 

'da f^Y e -^H= f rf S / ,/; " 







du J L V (is 



where 

Jds = e (3u ^du and r{£/(i/)} = / e^^du (A5) 

Jo 

The new final time t{L / "( , u)} is not fixed anymore, but depends on the realization of 
random potential U(u). To take into account this constraint, we can insert the identity 

1 = jf°° dr s(t- J* eW^duj (A6) 

to obtain 

1 [\ ((3dU „ \ 2 

/•OO nx{t)=X r — — — / I V 2D [I 

X t (x)= dr Vx(u) VU(u)e 4D Jo V 2 cfo / 

io Jx(0)=0 Ju(o)=o 

1 /" T ( dx\ 2 
~4DJ dS W S (t- J*e? u Mdu^ (A.7) 

Let us now perform the Gaussian path-integral on a; (it) 

,x(t)=x --7T / (is 

%(«)e 40 -A) 

c(0) 

to get 



dec \ 2 x ? 



/ Vx(u)e AD Jo \dsj 1 4D t ( A8 ) 

Jx(o)=o v47rDr 



00 x2 

X t (x) = / dr— i e~4Dt ^(r) (A9) 

Jo VAtvDt 
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where 

1 /•* (pdU „ V 

Mt)= VU(u)e AD Jo \2 dv J sir- e^^du) (A10) 

Ju(o)=o \ Jo J 

is by definition the probability distribution of the functional r t = J * e^ u ^du. 

Eq (A. 9) is just the translation in terms of probability distributions of the identity in 

law between the process x(t) and a linear Brownian motion of stochastic time 7~t 



x(t) = V2D [ du 7](u) 
Jo 
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